Numerical approach for retention characteristics of double floating-gate memories 
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We report on a numerical investigation in which memory characteristics of double floating-gate 
(DFG) structure were compared to those of the conventional single floating-gate structure, including 
an interference effect between two cells. We found that the advantage of the DFG structure is its 
longer retention time and the disadvantage is its smaller threshold voltage shift. We also provide 
an analytical form of charging energy including the interference effect. 
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Floating-gate (FG) memories are widely used in com- 
puters because of their low cost and high density [l|, y] ■ 
FG memories have progressed rapidly through down- 
scaling using state of the art technology. However, sev- 
eral problems have been arising as a result of the progress 
of down-scaling of the FG structure to the nanoscale re- 
gion. In particular, the interference between FGs due 
to Coulomb interaction is emerging as one of the largest 
obstacles for FG memories [3-[a] ■ Stored charges in neigh- 
boring FGs interfere with one another, resulting in unde- 
sirable threshold voltage shifts in memory operations. In 
order to reduce this interference, complicated program- 
ming sequences are carried out in the current commercial 
FG arrays. 

In the case of locally charged materials such as dielec- 
tric materials, it is evident that the electric dipoles exist 
stably within mutual strong Coulomb interaction. This 
leads us to consider whether we can construct an "arti- 
ficial dipole" using a FG system. One of the candidates 
might be a stacked double floating-gate (DFG) structure 
Q , in which an additional FG exists between the FG and 
the control gate in the conventional FG array as shown 
in Fig. [TJa). Moreover, if the FG becomes as small as a 
quantum dot [3,[3j DFG can be used as a qubitQ, which 
is a basic element of a quantum computer [10|. There- 
fore, it is important to clarify the fundamental proper- 
ties of the DFG structure. The purpose of this paper is 
to numerically compare the retention time of the DFG 
structure with that of the conventional single FG (SFG) 
structure in Si/Si02 system, including an interference ef- 
fect between two cells. We clarify the unique transient 
behavior of DFG owing to the existence of the additional 
FG. In order to compare DFG with SFG impartially, 
we adopt the same equivalent oxide thickness (EOT) for 
both structures. We also compare read disturbs of both 
FG structures. Finally, we derive an analytical form of 
charging energy of DFG and SFG as a function of gate 
voltages and electron charges. 

Formulation. — We calculate transient behaviors of a 
cell where length L and width W of each FG and a 
distance between neighboring cells X-q are set equal as 
L — W = Xd = 23 nm and the height of all FGs is 
Z = 50 nm for two cases of oxide thickness (Fig. [T{b)). 
(We obtain similar results for the L = W =llnm 
case.) We take dielectric constants and effective mass 
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FIG. 1: (a) Double floating-gate (DFG) structure. Toxi, To^2 
and TcG are oxide thickness between lower FG and substrate, 
between two FGs, and between upper FG and control gate, 
respectively. qAi, qsi, etc. are stored charges in those capac- 
itances, (b) We consider two cases of device parameters. 



of Si and oxide Si02 as esi = 11.7, eox = 3.9, and 
TOgi — 0.19, TOox = 0.5, respectively. The barrier height 
of Si02 is $6 = 2.9eV. The capacitances are defined by 
Ca = esiLW/{TcG + lZ/2), Cb = esiLW/iTo^2 + jZ), 
Cc = esiLW/in^i + 7^/2), Cd = esiZW/^Xn + 
jL), Ce = esiZW/XE, Ch - esiZW/X^i and Ck = 
esiZW/XK with Xe - ^/(X^ 
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TcG are oxide thickness between lower FG and substrate, 
between two FGs, and between upper FG and control 
gate, respectively (Fig. 1). 7 = Cox/esi is a penetration 
effect of small FGs [iJl- Because we use Si02 for all tun- 
neling barriers, equal EOT means that the total thickness 
of barriers is the same, namely, T^^t^ - 

^SFG I T^SFG 

-f ox "T -f CG ■ , — , 

Transient calculation is carried out as follows [ll|. (i) 
For given initial charges Qoi > Q/3i (* = 1 and « = 2 indi- 
cate left cell and right cell, a and /3 indicate upper FG 
and lower FG, respectively), potential energies of FGs 
V^Qi, Vjii (i = 1,2) are obtained by solving the matrix 
equations: fll]. 
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FIG. 2: (a) Stored charge in "00" state (schematic of charge 
distribution of DFGs shown in inset) and (b) threshold voh- 
age shift AVth as a function of gate voltage for case I. (c) 
Program and erase characteristics of case I and (d) those of 
case II. 
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FIG. 3: Transient behavior of AVth (retention characteris- 
tics). (a)(b) AVth of (00) and single cell ("alone") states for 
case I ((a)) and case II ((b)). (c)(d) AKh of the left cell 
for (01) states, (i) Vba = 0; Vbi = -V'prg for DFG and 
Vgi = Vp.g for SFG. (ii) Vc^ = -Vp.g/2; Vg, = -Vp.g for 
both DFG and SFG. 



{i — 2,1 when i ~ 1,2). 14ub is a substrate bias 
and we set T^ub = 0. (ii) Once potential energies 
of FGs are determined, electric field applied on each 
oxide is calculated as the difference of potential en- 
ergies of FGs. For example, electric field between 
the stacked FGs is given by Eu = (Vg. — T4ub)/'?oxi. 
(iii) Current through each oxide is calculated by a 
direct tunneling model from applied electric field E 
as JiE) = AE'^exY>{-B[l-{l-ETo^/^b)^/'^]/E} 
with A = e^msi/(167r^?imox$fc) and B = 



4^/2^ 



ox$, 



3/2 



/{She) [1^. (iv) Then, new charge distribu- 



tion is obtained, namely, Qp- => Qp. + {Ju — J2i)dt with 
the current that flows through the lowest tunnehng oxide 
{Ju) and the middle tunneling oxide {J2i) during time 
dt. We repeat this calculation until charge distribution 
is stabilized by adjusting small time advance dt. To 
determine stored charge and WRITE/ERASE process 
for a given gate voltage Vpig' ^^ start from trial charge 
(10^^ C/cm^) and repeatedly apply Vprg and — Vprg a 
couple of times. The retention behavior is described 
under Vg = 0, starting from the stored charges. 

Numerical results. — First, we found that DFG is more 
stable when Tqxi = Tox2. For example, charge distri- 
bution of Toxi — Tox2 = 5 nm DFG begins to change 
later than that of Toxi = 6 nm and Tox2 = 4 nm DFG. 
Thus, we consider DFG with Tqxi = Tox2. This is be- 
cause charge distribution begins to change through the 
thinnest tunneling oxide. 

Figure[2](a) shows stored charges for programmed "00" 
states in case I. (We obtain a similar behavior for case 
II.) We define "0" state as a negative charge stored state 
(programmed state) such as Qa + Q/i < and "1" as a 
charge unstored state such as Qa + Qp > 0. For DFG, 
the upper FG stores negative charges similar to the FG 
of SFG, whereas the lower FG stores positive charges 
as if DFG constructs an "artificial electric dipole" . Fig- 
ure[2](b) shows the threshold voltage shift AVth for DFG 



and SFG. For given charges Qa and Qp, AVth is ob- 
tained by Eqs.(IT]) as a gate voltage shift when Ei^i = 0. 
For "00" state, AVth is given when Ei^i = Ei^2 = 
and, for "01" state, AVth is given only for the left cell 
when _Ei,i = = Vg2. We can see that the magnitude 
of AVth of the DFG is one-fourth smaller than that of 
the SFG. This is because positive charge and negative 
charge cancel electric fields with each other and the elec- 
tric field outside the dipole structure of DFG is weaker 
than that outside a single charge structure of SFG. This 
indicates that DFG is less appropriate for memories us- 
ing multi- levels than SFG. In Figs.[2l!^c)(d), we show that 
the WRITE/ERASE speeds are almost the same in both 
structures. This is because the speed is mainly deter- 
mined by the same Tqxi . The peaks of DFG in the figures 
originate from different changes of Qa and Qp and ap- 
pear when the sign of Qa + Qp is changed. Because the 
charge distribution in our model is symmetric for Vg > 
and Vg < 0, the smaller memory window of DFG in 
Figs.[2]Jc)(d) corresponds to twice the AVth in Fig. [2)^b). 
Figures [3] (a) (b) show retention characteristics, that 
is, transient degradations of AVth of DFG and SFG for 
"00" and "alone" states. We can see that AVth of DFG 
exceeds that of SFG at ~ 10^° sec for case I and at 
^ 10^ sec for case II. Thus, retention time of DFG 
in "00" state is longer than that of SFG. The peak of 
DFG AVth appears when Qa + Qp begins to decrease. 
Note that AVth of "alone" state starts from negative re- 
gion. Thus, the interference between cells is effective 
for DFG. Programming voltage for (01) state is differ- 
ent from that for (00) state. In Fig. [3] (c)(d), negative 
voltage is applied to store negative charge for DFG and 
SFG of Vg2 7^ 0. These are results of the comphcated 
electromagnetic fields produced by Coulomb interactions 
among FGs. In any case, AVthS of "01" state in DFG 
become larger than those of SFG even in this thin Tqxi 
(< Snm) region. (We have similar relation between DFG 
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FIG. 4: Read disturb at Kcad = 4 V starting from "00" states 
for case I. Koad is applied to the left cell. Note that AVthS 
for SFG cells do not change during this time scale. 



The charging energy of the system can be expressed by 
summing charging energy of all capacitors such as C/ch = 
'^iqf/i'^Ci) — J2i' Qi'^v where qi shows charge of each 
capacitor and qii shows charge of capacitor that connects 
to gate voltage (see Fig. 1). By using Lagrange multipli- 
ers similar to Ref . [9| , we have the charging energy of two 

DFGs as C/dfg = Ui + [/„ - ELiI^a.^J, + Ch, V(?,}/2, 
where 



DIWI^ + D2'WI^ + 2CyWg^W, 



Ui = 



UU = {Ca.vl + Ca.vl + 2CD«g,«gJ/(2A) (2) 



2[DiD2 - C2] 

"gi ' '^u-i"g2 ' "^D'^Si ''32; 
-iDFG 



with vg^ = Ca,Vg^ + Ch,Vg. - Qr and 



and SFG for T^J^ = TD^g 



and TSfg 



■^0x2 



= T^/G = 7nm, Tgi^ = 13 
20nm [Ij].) Once dipole is formed in DFG, 
its surrounding electric field is weaker than that of SFG. 
Thus, DFG weakly couples with its environment, result- 
ing in longer retention time. The interference between 
two cells is considered to stabilize the charge redistribu- 
tion of horizontal directions. 

Read disturb. — Figure |4] shows retention characteris- 
tics in which gate bias Vqi = VJ-ead = 4 V is applied on 
the left cell while Vqs = 0, assuming that the left cell is 
read by V^oad (read disturb process). We obtain a similar 
behavior for case II. As can be seen, DFG changes earlier 
than SFG, but because the reading process is carried out 
in millisecond order, this weakness does not affect the 
practical usage of DFG. 

Charging energy. — In the near future of the down- 
scaling of FGs, the number of electrons in the FGs is 
reduced and countable, resulting in the region of single- 
electronics 15|. In this region, an analytical form of 



charging energy as functions of electrical charges and 
gate biases is required to describe an experimental sta- 
bility diagram of electron charge distribution as shown 
in Ref. Jil . Here we provide an analytical form of the 
charging energy of two cells (Fig. [T]) by using a capaci- 
tance network model. 



+ {[CdCb, +Ca,CEj]/AK, -Ql^fG, (3) 
including C^, = Ca. + Cb. + Cd + Ce. + Cr,, Cb, = 

Cb; + Cci + Cg + Ce-- + Ckj, a = CaiCa2 — Cq, 

A = a, - [Ca.Cl + Ca.Ci, + 2CdCe,CbJ/A, and 

Cy = Cq + [CaiCBaCEa + C'a2C'BiC'Ei + CbjCbsC'd + 

CeiC'esC'dI/A. Note that Usfg has the same form as 
Uu with capacitances replaced by those of SFG such as 

Wg^ and 



Cl„ C%i etc. For "00" state {wg 
Ugj), we have. 



ugi 



Ui = wy{D - Cy), Uu = vlliCa - Cd) (4) 

For SFG, we have C/sfg = fg/(Cf - Cg) - (Ca + Ch)1^J 
where Cl ^ C\-V C% + C^ + C\ + Cf and Vg = {C\ -K 

cDyG-QT""- ^ ' _ 

In summary, we numerically showed that interfering 
DFGs have a longer retention time than SFGs, although 
DFGs have the disadvantage of smaller threshold volt- 
age shift and smaller memory window. Owing to the 
existence of additional FG, DFG shows unique transient 
characteristics, forming an artificial electric dipole. 

The authors thank A. Nishiyama, J. Koga, S. Fujita, 
N. Yasuda and A. Kinoshita for useful discussions. 
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